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Abstract 

We study non-Gaussianity, the spectral index of primordial scalar fluctuations 
and tensor modes in models where fluctuations from the inflaton and the curvaton 
can both contribute to the present cosmic density fluctuations. Even though simple 
single-field inflation models generate only tiny non-Gaussianity, if we consider such a 
mixed scenario, large non-Gaussianity can be produced. Furthermore, we study the 
inflationary parameters such as the spectral index and the tensor-to-scalar ratio in 
this kind of models and discuss in what cases models predict the spectral index and 
tensor modes allowed by the current data while generating large non-Gaussianity, 
which may have many implications for model-buildings of the inflationary universe. 



1 Introduction 



It is widely believed that cosmic microwave background anisotropics and large scale struc- 
ture of the universe we observe today originate from fluctuations generated during the 
time of inflation. It is often assumed that the quantum fluctuation of the inflaton is re- 
sponsible for that. During inflation, the inflaton slowly rolls down its potential, which 
gives almost scale-invariant primordial fluctuations and some gravity waves can also be 
generated. Usually the primordial fluctuations are characterized by the scalar spectral 
index and the tensor-to-scalar ratio. Once the potential for the inflaton is given, one can 
predict these quantities to compare with cosmological observations. Current cosmologi- 
cal data are very precise to severely constrain models of inflation and some models are 
considered to have been already excluded [H [2] . 

In fact, the considerations of the spectral index and tensor modes may not be enough 
to probe primordial fluctuations. The non-Gaussianity of the fluctuations can also be 
measured by cosmological observations and be used to test the scenario generating pri- 
mordial fluctuations. In particular, there has recently been reported that non-Gaussianity 
is detected in the cosmic microwave background almost at 3a level [3]. Since simple 
single-field inflation models predict very small non-Gaussianity, if the evidence for large 
non-Gaussianity is mounting in the future, these models which at present are judged to 
be successful as regards the amplitude and scale dependence for primordial curvature 
fluctuations may suffer from a great difficulty. 

However, notice that another source of fluctuations other than the inflaton can con- 
tribute to the cosmic density fluctuations. Among such possibilities, the curvaton mecha- 
nism [H El E] and modulated reheating scenarios [71 E], for example, have been proposed 
and some of their observational consequences including issues of the non-Gaussianity have 
also been investigated. In particular, it is discussed that large non-Gaussianity can be 
produced in these scenarios. It should be reminded that cosmic density fluctuations can 
originate from one of these mechanisms including the inflaton or from some combination 
of these. Thus, even though almost perfect Gaussian fluctuations are generated in simple 
inflation models where they originate only from the inflaton, large non-Gaussianity can be 
produced for total by adding the contribution from another generation mechanism. One 
of the purposes of this paper is to study the issue of non-Gaussianity in models of this 
kind, paying particular attention to a mixed scenario where fluctuations of the inflaton 
and the curvaton can both contribute to primordial fluctuations. In investigating the non- 
Gaussianity, we focus not only on the bispectrum but also on the trispectrum which is a 
potentially useful observable in the near future. 

It should also be noticed that the predictions of the amplitude, the spectral index 
and tensor modes can be affected when another source of fluctuations is introduced in 
the model. Interestingly, some models of inflation which are disfavored by the data can 
be liberated by assuming the curvaton mechanism to work because of such modifications 
il[ini[IIl[ni[I31[Il[ISl[IS]. In particular, in Ref. [B], it has been studied in what 
cases models of inflation can be relaxed by adding fluctuations from the curvaton in some 
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detail assuming some concrete inflation models focusing on the spectral index and the 
tensor-to-scalar ratio. In any case, since we can expect much more precise measurements 
of these quantities in the future experiment such as Planck, we should also investigate 
the predictions for the scale dependence of primordial curvature fluctuations and tensor 
modes in the mixed scenario as well as that for non-Gaussianity. Another aim of this 
paper is to analyze this issue extending the work of Ref. [15] to include the case where 
the curvaton is always subdominant component in the history of the universe, which was 
not considered in Ref. [15]. Then we study the effects of the curvaton on inflationary 
parameters such as the spectral index and the tensor-to-scalar ratio and compare with 
recent observations of WMAP. Furthermore, we also discuss in what cases/models large 
non-Gaussianity can be generated satisfying the constraints on the scale dependence and 
tensor modes of primordial fluctuations in the scenario. 

The structure of this paper is as follows. In the next section, we discuss the formalism 
to study observational quantities in models with mixed fluctuations from the inflaton and 
the curvaton. We derive the expressions for the scalar spectral index, its running, the 
tensor-to-scalar ratio and the non-linearity parameters for such mixed models. Then, in 
Section [3|, we investigate those inflationary parameters assuming several concrete inflation 
models. First of all, we study the scalar spectral index and the tensor-to-scalar ratio 
for each model and then discuss whether they are compatible with current observations. 
For inflation models which are considered to have been already excluded by the data, 
we discuss in what case the curvaton can liberate the model. Specifically, we discuss in 
what cases non-Gaussianity can be large. The final section is devoted to conclusion and 
summary of this paper. 

2 Formalism 

2.1 6N formalism and some definitions 

In the following, we adopt the 6N formalism [TTJ [181 UHl 120] to calculate the primor- 
dial power spectrum and some non-linearity parameters. Thus we briefly review the 6N 
formalism here. 

In the 6N formalism, the primordial curvature perturbation ( on the uniform energy 
density hypersurface at the time t = tf is equal to the perturbation in the local expansion 
defined with respect to an initial spatially flat hypersurface. On sufficiently large scales 
where the spatial gradient can be neglected, the local expansion is well approximated by 
the expansion of the unperturbed universe. 



where H is the local Hubble expansion and is some time during inflation. Then the 
primordial curvature perturbation can be expressed as 




(1) 




(2) 
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where is the expansion in the background spacetime: 

N= r H{t)dt. (3) 



If we take as a time when the cosmological scale crossed the horizon scale during 
inflation, then N{tf,t^,x) becomes a function of the scalar field at the time of horizon 
crossing. Hence ( can be written as, 

c(t/) = nm: + iNabSrMt + iN^bcsrM'Mt + ■ • • , (4) 

where the summation is implied over repeated indices which label the scalar field. Here 
50" is the perturbation of the scalar field 0" on the fiat slicing at the time of horizon 
crossing. In the following, when the asterisk * is shown in the subscript, it indicates that 
the quantities are evaluated at the time of horizon crossing. Na, Nab, and Nabc are given 
by 

" ~ or' " ordct)^' " drd(t)'d(t)-' 

For the purpose of this paper, we include the terms up to cubic order in the perturbations 
of the scalar field. If we choose tf well after the reheating, then C,{tf) gives the primordial 
adiabatic perturbations. 

Now we write down the expression for the primordial power spectrum which is defined 

as 

(C,^^Cg,) = (2vr)Vc(A;05(^i + 4). (6) 
By using Eq. (jl]), we can express with the fluctuations of scalar fields as 

P^{k) = NaN,P''\k), (7) 

where P°-^{k) is the power spectrum of scalar fields: 

= {27rf P^\k,)6{h + h). (8) 

In this paper, we consider two scalar fields, the infiaton and the curvaton which are denoted 
as and a and their fluctuations are indicated as 50 and 6a respectively. Thus the indices 
for scalar fields are understood to represent either of these fields in the following. The 
fluctuations at the time of horizon crossing 50* and 50"* are assumed to be uncorrelated 
random fields with the same amplitude. Hence we have 

P'''{k) = P{k)6'''='^(—] 5'^^ (9) 

k'^ \2ti J 

Then at leading order in fluctuations of the scalar fields, the power spectrum can be written 

as 



n(fc) = ^n(fc) = A^.Ar»— (^) . (10) 
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The expression for Na will be discussed in the next subsection. 

Now we discuss quantities which represent non-Gaussianity, the bispectrum and trispec- 
trum. The definition of the bispectrum is given by 

iCkSk^CkJ = (27r)^i?c(/ci, h)6{h + h + h). (11) 
As is the case with the power spectrum, the leading order bispectrum can be written as 

B^{h,h,h) = NaN,N,B-''{k,,h,h) 

+NM,Nd {P''\h)P'\k2) + P''\k2)P'\h) + P''\k^)P'>\k,)) , (12) 

where B"-^"^ is the bispectrum of the scalar fields 

= {2Txf B'^'^k,, k2, ks)6{h + h + h). (13) 

The trispectrum is defined as 

i(k,Ck,(k,Ck,) = W'Tdh, k2, k,, k,)S{h + h + h + h), (14) 

and can be given with the power spectrum P"'', the bispectrum 5°'"^ and the trispectrum 
T"-'^'^'^ for the scalar fields, 

T^iki, k2, ks, k,) = NamN,NdT''''"'{ku k2, ^3, k^) 

+NatN,NdN, {P-'{k{)B'''\k^2, k^) + ll perms.) 
+NabN,dNeN; {P'\k^^)P^\k;)P''f {k^) + 11 perms.) 
+NabcNdNeNf {P''\k2)P'\k:i)P'f{k^) + 3 perms.) , (15) 

where kij = \ki + kj\ and T"*'^'^ is given by 

{6<P15<P15,P16<PI) = {2'nfT'^"^\k,, ^2, ^3, ^4)5(^1 + k2 + h + h)- (16) 

Expressions for Babe and Tatcd were provided in [21] and [22] for general multiple-fields 
slow-roll infiation (see also [23]). In the following, we will neglect Babe and Tabed because 
the inclusion of these terms only gives corrections of slow-roll order to the primordial 
bispectrum and trispectrum, which are far below the observational sensitivity expected 
for the Planck satellite. Then B^ and can be written as 

Bdku k2, k,) = ^/nl (Pc(A;i)n(^2) + Pdk2)Pdk3) + Pciks)Pc{k,)) , (17) 

T(;{ki, k2, ks, ki) = tnl (Pc(^13)^c(^3)^c(^4) + H perms.) 
54 

+7^9nl iPdk2)Pdh)Pdh) + 3 perms.) , (18) 



4 



where /tvl, tnl and qnl are constant parameters which are given in terms of the deriva- 
tives of the number of e-folding with respect to the scalar fields as [25| [26] , 



6 _ KmN '^" 
5 (iVciV' 



/nl = (19) 



These expressions are valid for the case where Babe and Tabcd are negligibly small. Fur- 
thermore, an interesting inequality can be obtained with the help of the Cauchy-Schwartz 
inequality [27] : 

rNL > |/^L- (22) 

This inequality holds for the scenarios such as the curvaton and the modulated reheating 
scenarios where the leading non-Gaussianity comes from super-horizon evolution. 

To calculate the primordial power spectrum and the non-linearity parameters, we need 
to know the number of e-folding as a function of the field value of the scalar fields, the 
inflaton and the curvaton, from the time when the cosmological scales crossed the hori- 
zon to the time when the universe becomes radiation-dominated connected to big bang 



nucleosynthesis (BBN) 1^. In the next subsection, we discuss this issue in detail. 



2.2 Background dynamics and the number of e-folding 

Here we give an explicit expression for the number of e-folding in models with the inflaton 
and the curvaton. For the potential of the curvaton, we take a quadratic potential, 

U{cr) = \mla\ (23) 

where mo- is the curvaton mass. We assume that the energy density of the curvaton 
is subdominant during inflation and at least until the time when the inflaton decays 
into radiation. In addition, the the curvaton mass is assumed as mg- <^ i/jnf with ifjnf 
being the Hubble parameter during inflation. In this case, the curvaton field almost stays 
at the initial value during inflation. After inflation ends, the inflaton oscillates around 
the minimum of its potential for some time, then decays into radiation. The curvaton 
also begins to oscillate at around a time when H = m^j and decays into radiation. In 
the following, we assume that the curvaton decays long after the onset of the curvaton 

Since we consider the scenario with the curvaton, the universe could have experienced the dominance 
of radiation twice after inflation ends: radiations from the decay of inflaton and curvaton. Radiation- 
dominated epoch here is meant to be the radiation-dominated phase after the curvaton decay. 
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oscillations, which is equivalent to assume that T„/m„ <^ 1 with T„ being the decay rate 
of the curvaton. 

Here we give a rough sketch of the thermal history of the universe in this scenario. In 
fact, depending on the initial amplitude of the curvaton field, the thermal history after the 
inflaton decay can become different. When the initial value of the curvaton (Xin is small 
enough, typically as o-jn ^ A/pi, the curvaton begins to oscillate around the minimum 
of the potential during radiation-dominated epoch due to the inflaton decay. With the 
potential Eq. fl23|) . the energy density of the curvaton decreases as p^- oc when a 
oscillates. On the other hand, the energy density of radiation decreases as prad a~^, 
thus the curvaton can become dominant if the curvaton oscillates around the potential 
minimum long enough before decaying into radiation. If the decay rate of the curvaton is 
not so small, it can decay before it becomes a dominant component. In any case, when the 
rate of the Hubble expansion H becomes comparable to the decay rate of the curvaton, 
the curvaton decays into radiation. After the curvaton decays into radiation, the universe 
is radiation-dominated and is led to BBN epoch. 

When the initial value for the curvaton field is large enough, typically as (Jin ^ Mpi, 
the curvaton almost stays at the initial value and drives the second inflation even after the 
onset of the radiation-dominated epoch due to the inflaton decay. The thermal history 
in this case is as follows: after the inflation driven by the inflaton ends, the inflaton 
decays into radiation, then the universe becomes radiation dominated. When the initial 
amplitude of the curvaton is large, a stays at the initial value until long after the inflaton 
decays into radiation. Since the energy density of such a scalar field is constant, at some 
epoch, the energy density of the curvaton dominates to give the second inflationary epoch. 
This situation is similar to that of a double inflation model [28]. After the curvaton drives 
the second inflation, it begins to oscillate around the minimum of the potential when the 
Hubble parameter decreases as if ~ mg.. Then the curvaton decays into radiation when 

Since we would like to calculate the primordial power spectrum and some non-linearity 
parameters during the radiation-dominated epoch after the curvaton decay, we need to 
know the e-folding number from the time when the cosmological scale exits the horizon 
during inflation to the time after the decay of the curvaton. We denote the number of 
e-folding during these epochs as A^tot- For later convenience, we divide Atot into several 
parts as 

A^tot = A^inf + A^d + A^R (24) 

where Ainf represents the number of e-folding during inflation driven by the inflaton. A^d 
is for the epoch from the end of inflation to some time after the decay of the inflaton. A^r 
refers to the e-folding number from that time to the time well after the curvaton decay. 
Ar is relevant to the curvaton dynamics. 

Now we make detailed discussion on A^inf , A^d and A^r in order to calculate the primor- 
dial power spectrum and the non-linearity parameters in the 5N formalism. We do not 



6 



discuss Nd here since it is irrelevant to the primordial fluctuation'^^. However, when we 
relate the cosmological scales observed today to the time of horizon exit during inflation, 
we also need Nd, which we will discuss in the subsection 12.61 

First we start with iVinf, the number of e- folding during inflation. It is determined by 
the dynamics of the inflaton whose equation of motion is 

4> + 3H<P + V^ = 0, (25) 

where a dot represents the derivative with respect to the cosmic time and = dV/dcf). 
During the inflation, the so-called slow-roll (SR) approximation is valid. Hence A^inf can 
be written as, by using (f) ~ — V^/(3if), 

TVi^f = / Hdt^ -—^ / —d<P (26) 

where t^, and te respectively represent the times when a scale exits the horizon and when 
the inflation ends. Similarly, (p^ and 0e indicate the scalar fleld values at corresponding 
epochs. To give a more concrete expression for A^inf, we need to specify the potential for 
the inflaton. 

Next we consider A^r. After the inflaton decays into the radiation, the universe is 
composed of radiation and the curvaton. Hence the background equations are given by 

Pr + AHpr = T,p,, (27) 
a + (3if + T„)d + mla = 0, (28) 

^' = T7lL(Pr + Pa), (29) 



3M, 



where and p^. = (1/2)(t^ + U{(j) are energy densities of radiation and the curvaton, 
respectively. To obtain N-^., we solve the above equations from the time t = to which 
corresponds to some time after the reheating epoch due to the decay of the inflaton with 
the initial conditions 

2 

Pr{to) = Pro, Cr(to) = CT*, Cr(to) ^ (30) 

6H [to) 

to the time tf when the Hubble parameter becomes much smaller than T^, i.e.,H(tf) <^ T^. 
Here we denote the initial amplitude for the curvaton fleld with an asterisk as a*. Remind 
that we are using an asterisk to represent that a quantity is evaluated at the time of 
horizon exit. Since a almost stays at the initial position during inflation, the initial value 



'^'^ In the case that the curvaton energy density becomes dominant during the regime of the inflaton 
osciUations, N^i depends on cr* so that it is also relevant to the primordial fluctuations. However, as stated 
in the paragraph below Eq. (|23p. we do not consider such a case in this paper. 
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for a here is the same as that at the horizon exit during inflation. By integrating Eq. ([2 
we can write N-^ formally as 

Nn = ]\og^ + ] log (1 + m., r,)) , (31) 

where F(cr*, mo-, To-) is defined as 



oo 



F(a„m.,r.)=/ c^ATe^^-^^, (32) 
Jo H{N) Pro 

and pf = Pr(tf). In general, we need a numerical calculation to evaluate F(cr*, mg.. To-). 
However, for some cases, we can integrate analytically the right hand side (RHS) of 
Eq. (132|) . which will be explicitly given in the next subsection. It should be noted here 
that, when the curvaton does not drive the secondary inflation and begins to oscillate dur- 
ing the radiation-dominated epoch caused by the infiaton decay, the value of F depends 
only on the combination 

P= / 33 

When the curvaton is oscillating, it can be regarded as a matter fluid as long as the 
background evolution is concerned. Then the universe can be regarded as a two component 
system which consists of matter and radiation. In this case, it has been shown that the 
solution for the system of Eqs. (!27|) . (!28l) and (!29l) depends only on the single parameter 
p defined above [29]. Thus F can be determined once p is fixed for such a case. 

When the second inflation driven by the curvaton occurs, the curvaton cannot be 
regarded as matter. However, even in this case, the value of F can be determined once a* 
and the combination r„/m„ are given. During the second inflationary phase, the slow-roll 
approximation is valid for the curvaton field as in the case during inflation. In this case, 
TVr can be written as 

where cxend is the value of a at the end of the second inflation. C represents the e-folding 
number from the end of the second inflation to the time well after the curvaton decay and 
does not depend on a^,. Thus, even in this case, once the combination T^/rricj is given along 
with cr*, the number of e-folding iV^ is determined. Hence, we define the combination of 
the variables as 

s ^ I^, (35) 

which will be used frequently instead of giving Fg- and rricr separately. As stated before, 
we assume s <C 1 in this paper. 

Having described the necessary formulae in the 6N formalism, we discuss primordial 
curvature fluctuations, tensor modes and non-Gaussianity in the following subsections. 
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2.3 Primordial power spectrum 

Now we discuss primordial power spectrum in models with mixed fluctuations from the 
inflaton and the curvaton. By using the 6N formalism, the primordial curvature pertur- 
bation ( can be calculated as 

where Q represents the a* dependent part of which is defined as 

Q^llog(i + F). (37) 

In actual calculations, we numerically obtain the function F in all cases. The primordial 
power spectrum for a model with mixed inflaton and curvaton fluctuations is given as 

1 fHV 



where Qa = dQ/da and we used Eq. Qj. The spectral index and its running are given by 
dhgVc 



k=aH 



dlogk 

1 + 2eM^,Ql 



4e - 277 

26- ^ , , (39) 



dus 



dlogk 



(40) 



'^2^i[v) ' ^"^P'T^' ^"^K"^- (41) 

Here the third parameter should be considered to be the second order in the slow-roll. 
During inflation, tensor modes are also generated as 



where e, 7/ and ^ are the slow-roll parameters which are defined as 

/ \ 2 



8 fH 
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When one extracts the information on tensor modes, the tensor-to-scalar ratio is usually 
used, which is defined and given in the case considered here as 

Vt 16e 



Vc l + 2eM2gr 
9 



(43) 



In the pure curvaton limit, 2eMpYQl. ^ 1, where the primordial fluctuations are com- 
pletely dominated by the curvaton fluctuations, ns, nrun and r reduce to 

- 1 = -2e, (44) 
rirun = -4e(2e-r/), (45) 

^ ^ 4qi' '''' 

We flnd that the spectrum of the scalar fluctuations becomes red-tilted, the running of 
the spectral index becomes independent of and the tensor-to-scalar ratio becomes 
negligible. 

2.4 Non-linearity parameters 



In the subsection 12. we gave the deflnitions of non-linearity parameters of /nl, tnl and 
(7nl in the 6N formalism. Here we give concrete expressions for them in models with mixed 
fluctuations from the inflaton and the curvaton. 

To obtain the bispectrum and the trispectrum, we need to expand 6N up to cubic 
order in 6(f) and 6a, which can be written as 



+Qa6a, + + \QaaJ(yl- (47) 

2 

Then, using Eqs. (fT9ll . ( I20l) and ( |2T]) . the non-linearity parameters are given by 

tJ^- = (1 + AiiQiT - " + ^^''^^'Q^Q'-] • (48) 

54 1 

25^^^ = (l + 26Afp^,g^)3 ["2'^ - + + 86=^M5Q^g_] . (50) 
In the pure curvaton limit, 2eMp[(5^ ^ 1, these expressions reduce to 

^/nl = (51) 

T-NL = = ■^/nl) (52) 
= (53) 
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10* 10"^ 10"'' 10""'' 10"* 10"' 10"- 10"' 1 lO' 10"* 10"' 10"'' 10"' lO"* 10"-' 10"- 10"' 1 lO' 

a, a. 



Figure 1: Plots of Q (top left), Q^, (top right), Qo-o- (bottom right) and Qo-o-o- (bottom left) 
as functions of cr* for the cases with s = 10~^ (red solid line), s = 10~^ (blue dotted line) 
and s = 10~^^ (green dashed line). Notice that Qo-o- and Qaaa can be negative. When the 
functions take a negative value, the absolute value is drawn with a thin line, cx^, is shown 
in units of Mpi. 



Once we pick up a model of inflation and calculate the values of Qcr, Qaa and Qaacn we 
can give concrete values for these non-linearity parameters along with Ug and r. Although 
an analytic expression can be obtained for some limiting values of p and a*, which is 
going to be discussed in the next subsection, a numerical calculation is needed to evaluate 
Qa, Qua and Qaua in general. In Fig. [H we show them as functions of cr* for several values 
of s. From the figure, we can obtain some idea of when the contribution from the curvaton 
significantly affects the primordial power spectrum and non-Gaussianity. 

2.5 Analytic expressions in limiting cases 

As previously mentioned, the value of F can be obtained by numerical calculations in 
general. However, we can express it analytically in some limiting cases. In the following, we 
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begin with the case where the initial amphtude of the curvaton a* is very small (a* <^ Mpi). 
In fact, this case can be further divided into two situations. The first case is that energy 
density of the curvaton once dominates the universe before the decay of the curvaton. 
The other case is that it can be always neglected compared to radiation during the whole 
history of the universe. These two cases correspond to p ^ 1 and p <^ 1, respectively. 
After discussing these cases, we consider the situation where the initial amplitude for the 
curvaton is large enough to drive the second inflation. Now we look at these cases in order. 



2.5.1 The case with <^ Mpi and p^ 1 

Here we adopt the so-called sudden decay approximation, which can give a good description 
for the case with p ^ 1. In this approximation, the curvaton decays suddenly when H = Ta 
at TV = Nfi- That is, we shall regard energy densities pr and pa to behave like the step 
functions at the time of the curvaton decay. All p^. is transfered to pr- Hence p^ increases 
by Po- at the time of the curvaton decay. Then, the equation for radiation energy density 
is written as 

'^+Ap, = p^{N,^)6{N-N,), (54) 

where the equation is written in terms of the derivative with respect to the number of 
e-folding and A^^. represents an e-folding number just before the time when H = T(j. The 
coefficient on the RHS is chosen so that the total energy density is conserved through the 
curvaton decay. Then, the equation corresponding to Eq. fl32|) can be written as 

Now let us introduce a new quantity q defined as 

a = (56) 
^ 4p,(A^d-) + 3p.(A^d-)' ^ ' 

which roughly represents the fraction of energy density of the curvaton to the total one at 
the time of the curvaton decay. With this parameter q, F((T*, mg-. To-) can be written as 

4 q 

F(cr*,m^,r^) = . (57) 

3 1 — g 

Since we can relate g to a,,,, we can obtain the analytic expression of Qo- as a function of 
a^,m^ and T^. For this purpose, let us write down the relation between q and a^, in a 
following way. The Friedman equation at the time of the curvaton decay is given by 

H\N^) = = ^ (pr(A^d-) + P.(A^d-)) . (58) 
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Denoting the e-folding number at the time H = rricr as Nm, this equation can be written 



as 



3M2 



(smlMy^ + p.(iV^)e-3^) , (59) 



where = — N„^. Furthermore, by using the relation 



3m2Mp2 pr{Nd^) 31 -g 

one can relate q and Pa{N^) as: 



_ 1 ( pANJ \ (3 + g)^/^ ^,3/, 
16 I mlM^^ i 31/2^2 ^- 



qf'\ (61) 



Now we write the energy density of the curvaton at = as Pa{N^) = {l/2)m'la'^al 
where a parameterizes the amplitude of the curvaton at the onset of oscillations relative 
to the initial one. The concrete value of a is given in Appendix |Al By using them, we 
obtain the relation between q and a^^ as 

"^= (3 + g)V4(i-,)3/4 - (62) 

By using Eqs. f l57|) and fl62l) . we obtain the expressions for the derivatives of F with respect 
to (T, which is necessary to evaluate the inflationary parameters, for this limiting 

Qa = Qa. = ^4(3 - 4g - 2g=^), Q.aa = 7^4(-18 + q + 2Qq^ + 6g=^). (63) 

In the limit = 1, which is equivalent to p ^ 1, Q„ becomes 

= (64) 

This agrees with the result in [13] where asymptotic form of Q„ for a* <^ Mpi is provided 
under the assumption that curvaton once dominates the universe before the curvaton 
decays. Hence the sudden decay approximation can well describe Q for the case p ^ 1. 
Using Eq. ( l63l) . the non-linearity parameters for p ^ 1 can be written as 

^ + 9^ , 

^?NL = — f-2er/ - + 2r^^ + ^6^^ ) . (67) 



-NL--^^^2.-,)V^e3^^ (66) 

1 + „2 



1 81 a' 

54 1 / . .0 . o 256 Ml 



9 o-J 



25^^^"" A , ,Ml\^ V ' ' ' ' ' 81 
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All the magnitudes of these non-linearity parameters become the largest when cr* ~ ^/eMpl. 
Substituting a* = ^/eMp\ to the equations above, we find that the maximum magnitudes 
of /nl, tnl and (7nl are all 0{1). Hence large non-Gaussianity is not produced for p ^ 1. 



2.5.2 The case with < Mpi and p < 1 

Now we consider the second case where the curvaton is always subdominant in the his- 
tory of the universe. Even in this case, the sudden decay approximation adopted above 
can give a good estimate to some extent, but its accuracy is 0(10%) compared to nu- 



merical calculation^^. However, another approach can give a better analytic expression 
for Q. Detailed descriptions for the derivation of the analytic formulae can be found in 
Appendix [HI Here we just give the expressions. 

Since here we are considering the case with p <^ 1, we expand F up to the second order 
in p, which is necessary for obtaining the formula for the third derivative of Q with respect 
to cr because the third derivative of p with respect to cr vanishes. After some calculations, 
we can show that the function F can be expanded as 

F{p) = l^^'P + Yii«V + 0{p% (68) 

where we neglected higher order terms in s assuming s <^ 1. Thus the derivatives of Q 
with respect to cr are given for this case as 

1 ~ 



= -^('T - 1)""^ (1 + OW) . (71) 



Here notice that, while Q„ and Qo-o- come from the leading order term in p, Q^aa comes 
from the next-to-leading order term in p. With the expressions above, we can explicitly 
write down the non-linearity parameters as: 



6 1 fn 2 e f ira^ p 



/nl = iJ^c^'^Df I ) (1 + 0{p)) , (72) 



1 / \ 2 / 4 



|*.^-(.-l)v/f"^^«»(^^^)(lH-0(p),, (74) 



Using the fitting formula, the sudden-decay approximation can give an accurate estimate for some 
parameter range |30] . 
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where the functions D f , and Dg are defined as 



^'W^IiT^- ^'(-)-(rf^- "'W^ofW' 

In deriving these expressions, we have neglected the non-Gaussianity coming from the 
infiaton fiuctuations because they are small quantities. In the pure curvaton limit, the 
non-linearity parameters given above become 

6 _ 24 1 288 1 54 288(7r - 1) 1 

r/nL - —7=—-, Tnl - 1^, T^fl'NL /— r „ • (76) 



Using Eqs. (!72|) . (1731) and (!74j) . we can derive a relation among the non-linearity pa- 
rameters specific to this case (i.e., p <^1). After a little arithmetic, we obtain 



48(7r-l) 

When p <^ 1, p can be written using /nl and tnl as 



vr 5184 ,_i//nl 



3 -l + 0{p). (77) 




Because of the inequality Eq. ([22]), the magnitude of p is smaller than that of f^l- Hence 
for the case of the very large non-Gaussianity (/nl ^ 1) which occurs only when p <^ 1, 
the RHS in Eq. fl77j) becomes very close to unity and the equation provides a simple 
consistency relation between the bispectrum and the trispectrum. Note that this relation 
holds not only for the pure curvaton model but also for the mixed model of the infiaton 
and the curvaton. Neglecting the second term on the RHS in Eq. fl77|) . we can further 
derive the inequality for /nl and (7nl as, assuming /nl being positive, 

< -M/,,. (79) 

OTT 

Thus if the very large non-Gaussianity is detected in the future, we can discriminate the 
mixed model of the infiaton and the curvaton from other scenarios that also generate large 
non-Gaussianity by using these consistency relations. 

Furthermore, notice that the argument x of the functions Dj, D^- and Dg corresponds to 
the ratio of the infiaton's fiuctuation to the curvaton's, i.e., x = 2eM^^Ql = Ccur/Cinf where 
Cinf and Ccur are the curvature fiuctuations generated from the infiaton and the curvaton 
respectively. Since all these functions take the maximum values at x ~ 1, non-linearity 
parameters become the largest when the fiuctuations from the curvaton are comparable 
to those generated from the infiaton with e and s being fixed. In this case, the initial 
amplitude of the curvaton becomes a* ~ (J*,ma,x ~ s/eMpi. If a^, is smaller than cr^^^max, 
the curvature perturbations are dominated by the infiaton fiuctuations which are highly 
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Gaussian. Hence when cr* is smaller than o"*,niax, the magnitudes of the non-linearity 
parameters take smaller values. On the other hand, if o"^, is larger than cr,,^inax; it becomes 
similar to the pure curvaton case, in which fluctuations from the curvaton dominate over 
those of the inflaton. In this case, /nl and (?nl are proportional to l/a1 and tnl are 
proportional to 1/ cr^. Hence, as we take the value of a* be larger than cr*,max, the magnitude 
of the non-linearity parameters becomes smaller. Thus if we vary a,, under the condition 
that e and s are fixed, in which we are changing the amplitude of Ccur with C\ni being fixed, 
the non-Gaussianity becomes the largest when cr* ~ o"^,max; i-e., when the contributions 
from the inflaton and the curvaton are comparable. On the other hand, if we instead vary 
e under the condition that s and a* are fixed, where we are changing the amplitude of 
Cinf with ^cur being fixed, the non-Gaussianity becomes the largest when the contributions 
from the inflaton are absent. 

The maximum values of /nl, tnl and the minimum value of (^nl for the variation of a* 
are given by 

6 1 /tF 2 e 1 

5/n,.,.„ = ^^2«^~— , (80) 

TTo;^ 1 

rNL,„„ = — - ~ ^. (81) 
54 4(71- 1) /7 2 e 1 



where Pmax is the value of p with a* = cT*^max, which is much smaller than 1. We find that 
the magnitudes of these non-linearity parameters have the same p dependence as those in 
the case of the pure curvaton, i.e., /nl and (7nl are roughly given by the inverse of the 
fraction of the curvaton energy density to the total one at the time of the curvaton decay 
and Tnl is roughly given by /nl 



#4 



As a final remark of this subsection, we comment on the relation tnl ~ /nl which 
ids for = a^^max- 1 
From Eq. (175|) . we have 



holds for (7=1, = (J^^max- In fact, this relation breaks down when a* deviates from a^^max- 



DJx) 1 , , 

T^H = 1 + -• 83 

Hence when the curvaton contribution to the curvature perturbations is subdominant, 
i.e., X <S 1, Tnl is enhanced by a factor of x~^ compared with /^l. This shows that 
the trispectrum as well as the bispectrum could also be important in detecting the non- 
Gaussianity of the primordial perturbations. As we will see in Section [3] where /nl, S'nl 
and Tnl are calculated for various inflation models in the mixed inflaton and curvaton 



In Ref. |31] . by requiring that the homogeneous energy density of the curvaton should be larger 
than its gradient energy density and the decay rate T„ is at least of order that given by the gravitational 
strength, an upper bound on /nl was provided as a function of the tensor-to-scalar ratio in the pure 
curvaton case. 
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scenario, there are parameter regions where /nljS'nl ^ C'(l) but tnl ^ 1. In such a case, 
the leading non-Gaussianity comes not from the bispectrum but from the trispectrum 
through the Tnl terms. 

2.5.3 The case with cr^ > Mpi 

When the initial amplitude is large enough, the curvaton field can drive the second inflation 
after the first inflation caused by the inflaton. In this case, TVr is calculated from Eq. ( IMl) 
and is given by 

N^^^y.-al) + C. (84) 
Thus, in this case, Q^-, Q^j^ and Qaaa become 

Qo" = 2 J^2 ' '^o-cr = 2j^^2 ' Q(raa = 0. (85) 



Then, the non-linearity parameters are given by 



1 .9 1 , 



rNL = — ^^|(2e-,)^ + -e3-^^ (87) 



^9nl = ^—-^{-2er] - + 2r]^). (88) 



25' 



2 ^ 



As for /nl and Tnl, the contributions from the curvaton become the largest when a* ~ 
Mp\/^/e. Substituting a* ~ Mpi/^/e to Eqs. ( l86l) and (IHTl) . we find that the maximum 
contributions of the curvaton to /nl and tnl are 0{e) and O(e^), respectively. As for g^i^, 
addition of the curvaton always reduces the magnitude of (7nl- Hence /nl, tnl and ^^nl 
are at most (9(e, 77), 0(e^, e?7, 77^) and 0{e'^,er],r]'^,^^), respectively. 

2.6 Corresponding Scales 

To compare the prediction for the primordial curvature fluctuations and non-Gaussianity 
with observations, we need to specify when the present cosmological scale exited the 
horizon during inflation. Since k^^ = a^H^ holds when the scale with the wave number 
crossed the horizon, the reference scale /cref where we probe the primordial fluctuations at 
the present time is related to that at the horizon crossing during inflation as 
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where ao and Hq are the scale factor and the Hubble parameter at present. The ratio of 
a^: to flo ill the RHS can be divided into several parts as, 

^rcf ^end (^<f>reh (90) 

Here Cend; fl^rch and af are the scale factors at the times when the inflation driven by the 
inflaton ends, the inflaton decays to reheat the universe, and some time after the curvaton 
decays into radiation, respectively. For the definiteness, we take H ~ 10~^ro- as the time 
corresponding to a/ so that Eq. ( 13T1) gives a good description. By taking logarithm of both 
sides and reminding that the number of e-folding during inflation is A^inf = log(a*/aend), 
A'inf can be written as 

iV,.f = - log ^ + log ^ + log ^ + log ^ + log ^. (91) 

For the first term, we take k^d = 0.002 Mpc~^ in the analysis. The second term in the 
RHS corresponds to — A^^^^ in Eq. ( IMl) . i.e., the e-folding number from the end of inflation 
to the decay of the inflaton to the radiation. When the potential of the inflaton near the 
minimum is written as V (x (j)", its energy density decreases as oc (2-6^/(0+2) |32j_ xhus 
the second term can be written as 

Oend a + 2 PSreh 01 + 2 

log = — — log ^— ~ — — log — ^, (92) 

flt/ireh Oa Pond -Hcnd 

where pend and p^reh are energy density of inflaton when inflation ends and that of radiation 
at reheating. For the second approximate equality, we use the sudden decay approximation 
for the reheating F,^ = H{t^reh) with being the decay rate of the inflaton. For the 
analysis in the subsequent section, we assume = 10~^ifend for definiteness. 

The third term corresponds to —N^ in Eq. flM|) . The effect of the curvaton on the 
total e-folding number is contained in this term (see Eq. fl3T]) ). 

As for the fourth term, we assume that no more entropy is produced after the curvaton 
decays. Thus this term can be rewritten by using the conservation of the entropy density 
per comoving volume. Since the entropy density is given by s = {27T/45)g^sT^ with g^,s 
being the total number of effective massless degrees of freedom, we have the following 
relation, 

log— = -log— = -log -3. (93) 

ao 3 Sf 3 g*sfTf 

For g^^sf at the time of aj, we take g^^sf = 100. 



3 Observational quantities in the mixed models with 
the inflaton and the curvaton 

Now we discuss primordial curvature fluctuations and non-Gaussianity in the mixed mod- 
els with the inflaton and the curvaton adopting some concrete inflation models. Once 
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the potential for the inflaton is specified, using the formahsm developed above, we can 
make firm predictions for the scalar spectral index, the tensor-to-scalar ratio and the non- 
Gaussianity for the model. The issues of the modification to the spectral index, its running 
and tensor modes have been studied in the case where the energy density of the curvaton 
should dominate the universe at late time [TU [151 CB] • Here we also include the case where 
the curvaton is always subdominant in the whole history of the universe. Furthermore, we 
give the predictions for the non-linearity parameters /nl, tnl and qjsil for infiation models 
in the scenario. In the following we consider chaotic infiation, new infiation, and hybrid 
infiation in order. 

3.1 Chaotic inflation 

10^ 



1 

10-^ 



10'' 



10''^ lO'" 10'** 10'' 10'" 

s 

Figure 2: Contours of the e-folding number N in the s-a* plane. Here we assumed the 
chaotic infiation model with n = 4 for concreteness. a* is shown in units of Mpi. 



Chaotic infiation [33] is described by the potential of the formj^ 

where A is a model parameter which can be fixed by requiring that the fiuctuations have 
the right amplitude to be consistent with observations. To fix the model parameters, we 
impose the WMAP normalization ^2] for the primordial curvature fiuctuations given in 
Eq. ( l38l) to fix the value of A. For different models of inflation discussed in the following, 
we also use the WMAP normalization to flx one of the parameters in the inflaton potential. 

T^^This type of a simple polynomial potential can be realized in supergravity [Ml[551[551[571[551[551HniHT] . 
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Figure 3: Contours of the ratio of the contribution from the curvaton fluctuation and the 
inflaton fluctuation Ccur/Ckif ^''^'^ shown in the s-a* plane. Here we assumed the chaotic 
inflation model with n = 4 for concreteness. For reference, we also show contours of p 
which enables us to see the region where an analytic expression is valid provided in the 
previous section. 



Here we consider some versions of chaotic inflation by taking several values of n which is 
assumed to be positive even number. With the potential Eq. ( l94l) . the slow-roll parameters 
in this model are given by 

The number of e-folding during inflation can be written as 

Nrni = ^^y.-<Pendf, (96) 

where 0end refers to the value of the inflaton at the end of inflation. In general, we 
can assume 0* 2> 0cnd- Thus we have a simple relation between 0* and iVjnf. When 
fluctuations of the inflaton alone are assumed to exist, the spectral index, its running and 
the tensor-to-scalar ratio are given by, using the number of e-folding A'inf, 

n,-l = -!1±^, (97) 



"■run oAr2 ' ^ 



n + 2 

Yi 

i^^ (99) 

^Vinf 
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Figure 4: la and 2(t allowed contours from WMAP3 alone are shown, which is generated 
from the chain provided in the webpage [l3]. 



When fluctuations from the inflaton alone are considered, the non-linearity parameters are 
very small, which are of the order of slow-roll parameters. An explicit expression can be 
obtained by setting the curvaton contributions such as and Qo-o- to zero in Eqs. ( l48l) . 
(gnD and dSOD. 

In the following, we consider the case with n = 2, 4 and 6, then discuss how the curvaton 
fluctuations change the predictions of the original inflation models. As we discussed in the 
previous section, the contribution from the curvaton fluctuations modifies Us^n^^.^ and r 
as Eqs. (HUI) and In addition, the curvaton affects the background evolution, 

which changes the number of e-folding during inflation. When the oscillating curvaton field 
dominates, its energy density decreases as pa oc a~^. If this phase lasts long, it reduces 
A^inf- When the curvaton drives the second inflation, the number of e-folding is drastically 
reduced as A'inf ~ 20-30. In this case, as can be read off from Eqs. ( l97|) . ( l98l) and ( l99ll . 
the spectral index becomes more red-tilted, its running becomes larger and the tensor-to- 
scalar ratio is also a bit larger due to the change of A^inf. Thus when one considers the 
effects of the curvaton, the modification to the background evolution is also important. 

To see how the curvaton parameters affect the e-folding number during inflation, con- 



tours of A^inf in the s-cr* plane are shown in Fig. For concreteness, we assumed n = 4 



For the range of s, we take 10"^^ < s < 10^'' in the analysis here. Although s is determined by 
the combination of the decay rate the curvaton mass as s = Ta/Tnu, to make our analysis general, we do 
not specify the values of Ta and rricr- However, one should keep in mind that there are some constraints 
on Ta coming from the bounds on the reheating temperature Treh- A possible upper bound for Troh 
comes from the consideration of the gravitino problem. To avoid the overproduction of gravitinos, the 
reheating temperature should be low enough as T^ch ^ 10^ GeV [331 US]. In addition, Treh ~ 10 MeV is a 
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to plot the figure. In the figure, a* is shown in units of Mpi. 

As is clearly seen from the figure, when the initial amplitude is large, A^inf is significantly 
reduced due to the second infiation driven by the curvaton. On the other hand, when 
(T* <C Mpi, dependence of iVinf on a^, and s becomes mild. The qualitative behavior of 
this result can be understood as follows. There are two possible sources that affect N-u^f 
due to the existence of the curvaton. The first one is that the curvaton adds extra e- 
folding number Q compared to when the curvaton is absent. As a result, A^inf must be 
replaced with Ai^f — Q. The second one is that if the curvaton fiuctuations contribute 
to the density fiuctuations, the energy scale of infiation must decrease so that the sum of 
the fiuctuations originating from both the infiaton and the curvaton satisfies the WMAP 
normalization. This effect is described by Q^. Hence A^i^f is a function of Q and Q^, i.e., 

iVinf = iVinf(Q,Qa)■ 
When 0"* ^ Mpi, Q depends only on (see Eq. (!84l) ). Hence the slope of the contour in 
Fig. [2] becomes zero. On the other hand, when a,, <^ Mpi, Q becomes a function of p while 
Qo- becomes a function of p and cr*/y^. Hence we can write A'inf as A'inf = A^mf (p, c"*/ 
When p <^ 1, the modification of Ajnf due to the existence of the curvaton mainly comes 
from Q„ which is a function of o'^/^/s at the leading order in p because Q is very small 
{Q = 0{p)). Hence for p <^ 1, the slope of the contour approaches 1/4. When p 1, 
both Q and can be important and the slope can vary along the contour. 

Furthermore, to see in what range of s and a* fiuctuations from the curvaton dominates 
over that from the infiaton, the contours of the ratio of the squared contributions from the 
curvaton and the infiaton, Ccur/Cif = '^^pi^Ql: are shown in the s-cr^ plane in Fig. [3l For 
concreteness, we again assumed the chaotic infiation model with n = 4. We also plotted 
contours of p = cr^/(Mp[-\/i), which enable us to see in what region the analytic estimate 
we discussed in the previous section can be valid. 

Now we discuss the infiationary parameters such as rig, r and the non-linearity param- 
eters in each case. First we study the case with n = 2. In this case, if we assume A^inf = 50 
and 60, the spectral index and the tensor-to-scalar ratio become {ns,r) = (0.96,0.16) and 
{ns,r) = (0.97,0.13), respectively. The running of the spectral index is negligible for this 
model. To compare these values with observations, we show 1 and 2a allowed regions 
from WMAP3 alone analysis in the Ug-r plane in Fig. |H By looking at the figure, we 
can see that the case with n = 2 is favored by WMAP3 data regardless of the curvaton 
contribution. 

If we include the curvaton, then the non-linearity parameters will become large and ns 



and r will also change Hence it is interesting to check that the inclusion of the curvaton 
does not spoil the successful values of rig and r, keeping the non-linearity parameters large. 

In Fig. [5|, contours of Ug and r are shown for the case with n = 2 when the curvaton 
is included in the s-cr* plane. With our assumptions explained in the previous section. 



possible lower bound not to spoil the success of BBN. (Detailed discussion on this issue can be found in 

[MIITIHHIIISIED].) 

Even when fluctuation of the curvaton is large enough to affect the values of these quantities, the 
running of the spectral index remains negligible in this model. 
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Figure 5: Contours of (left panel) and r (right panel) in the s-cx* plane for the chaotic 
inflation model with n = 2. In the absence of the curvaton contribution, the spectral index 
and the tensor-to-scalar ratio are Ug = 0.967 and r = 0.133 for A^inf = 60.3. 



the number of e-folding becomes as A^inf = 60.3, which yields the spectral index and the 
tensor-to-scalar ratio as (ns,r) = (0.967,0.133) without the curvaton contribution. We 
can see that, even when the curvaton is introduced to affect the primordial fluctuation, it 
is consistent with current observations. In Fig. El the non-linearity parameters are shown 
also in the s-cx* plane. When the curvaton drives the second inflation, the sign of /nl 
becomes negative and its size is 0{1) in the magnitude. When the curvaton dominates at 
least once in the history of the universe but there is no second inflation which corresponds 
to the region where a^, is small and s is large, the non-linearity parameters become positive 
but the size is also 0{1). As discussed in the previous section, non-Gaussianity becomes 
large when p <^1, which can be seen from the figure. This is why the long axis of contour 
ellipses lies along the line with p = al/M'^^'sf^ being constant. From the analysis here, it 
can be concluded that non-Gaussianity can be large when the curvaton also contributes 
to the primordial fluctuations without spoiling successful values of Ug and r. 

Next we discuss the case with n = 4. In fact, this model is on the verge of the 
exclusion even by the WMAP3 data alone. If we assume A^inf = 60, the inflationary 
parameters become (n^, r) = (0.95, 0.27), which are marginal at the current data. However, 
by introducing the curvaton in the model, the spectral index is modifled to be shifted to 
the scale-invariant one and the tensor-to-scalar ratio is suppressed. Thus it can liberate 
the model. In Fig. [71 contours of and r are shown in the s-a* plane. Here, in our setting, 
-^inf = 60.3, which yields (n^, r) = (0.950, 0.265) without the curvaton contribution. When 
the effects of the curvaton are small, the values of rig and r fall outside the allowed region. 
However, in particular, around the region where cr^, ^ Mpi and p <^ 1, the effects of 
the curvaton are signiflcant to make the inflationary parameters preferable values in some 
parameter space. To show in what cases the inflation model is liberated by the curvaton. 
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Figure 6: Contours of /nl (left panel), Tnl (center panel) and (^nl (right panel) in the s-cr* 
plane for the chaotic inflation model with n = 2. 



the excluded region from the WMAP3 data is represented by the shaded region in Fig. [H 
Namely, the region without the shade is allowed by the data. In the same figure, we also 
plot contours of the non-linearity parameters /nl,''nl and (^nl- Interestingly, there is a 
parameter space where the model can avoid the exclusion by the data while generating 
large non-Gaussianity. 

As for the case with n = 6, this model is completely ruled out by current observations 
without the curvaton. When A'^inf = 60 and 50, the spectral index and the tensor-to- 
scalar ratio are (ris,r) = (0.93,0.4) and (0.92,0.48). However, in the same way as the 
case with n = 4, this model can be made to be allowed with the help of the curvaton. 
In Fig. [HI contours of and r are shown. In our setting, A^i^f = 64.0, which yields 
(n5,r) = (0.937,0.375) without the curvaton contribution. In Fig. [TOl contours of the 
non-linearity parameters are shown along with the region excluded by the data in the 
s-(j^ plane. Since the original model of inflation is strongly disfavored by the data, the 
allowed region in Fig. [10] is smaller compared to the counterpart for the case with n = 4. 
However, there is a region where the model is relaxed to be allowed by observations and 
large non-Gaussianity can be generated. 

3.2 New inflation 

Among various possibilities of the potential for the new inflation model, we adopt the 
following form, which is motivated by the supersymmetric models with the discrete R 
symmetry of the group [511 [52], for definiteness: 

V = 
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Figure 7: Contours of (left panel) and r (right panel) in the s-cx* plane for the chaotic 
inflation model with n = 4. In the absence of the curvaton contribution, = 0.950 and 
r = 0.265 for iVinf = 60.3. 



During inflation, the last term in the RHS is irrelevant, thus neglecting this term, we can 
write down the slow-roll parameters for the potential above as 

r,--Mn-l)[^J-\^)\ (102) 

= 4n'(n - l)(n - 2) f^V'" " (^]' ■ (103) 
The number of e-folding is given approximately as 

2n(n-2) \MpJ \ v 

Since <^ 1 during inflation, the slow-roll parameter e is much smaller than rj in 

magnitude, which leads to the spectral index ri^ — 1 — 2?]. By using A^mf, this can be 
written as 

n — 1 1 , , 

- 1 ~ -2 -— . 105 

Furthermore, the running of the spectral index rirun — —2^"^ becomes 

n — 1 1 

'^-^ - "2-3^ T7^- (106) 



iVinf = . / f T^ ) (-) . (104) 
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Figure 8: Contours of /nl (left panel), tnl (center panel) and (^nl (right panel) in the 
s-(T* plane for the chaotic inflation model with n = 4. The shaded region is excluded by 
WMAP3 data. 



The tensor-to-scalar ratio r = 16e is very small since e is suppressed as 

.2 



In' 



V 



2 r 1 / \ n 2(n-l)/(n-2) 



_2n{n- 2)iVinf \Mpi 



(107) 



If we take n = 3 and 4, the e parameters are given by e = (f /Mpi)^/(72A'^^f) and e = 
4(w/Mpi)''/(277Vj^£), respectively. Thus even if we take v ~ Mpi, e is suppressed with 
some powers of A^inf which is usually 50-60. Hence both the tensor-to-scalar ratio and the 
running are very small in this model. The spectral indices for the cases with n = 3 and 
4 are = 0.93 and 0.95 when TVinf = 60. Thus this model is almost consistent with the 
current constraint from WMAP3 in terms of the spectral index and tensor mode. 

One may speculate that, adding a contribution from the curvaton, non-Gaussianity 
becomes significantly large like in the case of the chaotic inflation. However, this does not 
happen in the new inflation. This is because the slow-roll parameter e is very small in 
the new inflation models. Largest possible values for the non-linear parameters which are 
given in Eqs. fl72]) . fl75]) and flTll) indicate that these values are proportional to e. Thus 
inflation models with small e cannot generate large non-Gaussianity even by adding the 
curvaton. 

Furthermore, when e is very small, the effects of the curvaton on Ug and r are also 
small. As seen from Eqs. ( 1391) . (140!) and (l43l) . the effect of the curvaton on these quantities 
always appears in the combination of eQ^. Thus even when Qo- is large, this combination 
is small if e is negligibly small. To see this fact, contours of Ug for the case with n = 3 
and 4 are shown in Fig. [TTJ In both case, we take v = 10~^Mpi. Without the curvaton 
contribution, the e-folding numbers in our setting for n = 3 and 4 are given by A^inf = 50.3 
and 53.3 respectively. Then the spectral indices are = 0.920 and 0.944 respectively. As 
seen from the figure, the value of Ug does not change much over the whole parameter space 
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Figure 9: Contours of (left panel) and r (right panel) in the s-cx* plane for the chaotic 
inflation model with n = 6. In the absence of the curvaton contribution, the spectral index 
and the tensor-to-scalar ratio are Ug = 0.937 and r = 0.375 for A^inf = 64.0. 



in the figure. This indicates that the effect of the curvaton cannot be large for a broad 
range of the curvaton parameters. In fact, when the initial amplitude for the curvaton 
is large, Ug is modified a little because of the change of the number of e-folding due to 
the second inflation. However, in any case, the flgure clearly shows that the effect of the 
fluctuations of the curvaton is small for this kind of models. 



3.3 Hybrid inflation 



Hybrid inflation was originally proposed by Linde [53l and many types of the potential 
have been discussed in association with the gauge symmetry breaking [55l [561 EZ] • Here, 
we take the original types of the potential for simplicity. The implications of another fleld 
on other types of the potential like those in the supersymmetric grand unifled theory are 
discussed in the context of modular inflation [58] . 

We consider the following quadratic type of potential, 



V = a 



(V —a 



m r. q 
H d + 



Here we set the parameters to be w = 10~ Mpi, m = 2 x 10" Mpi, and g = 2 x 10" . The 
parameter a is chosen to flt the WMAP normalization. For > 0cnd = '^v/g = lOMpi, 
the trajectory with a = is stable and inflation takes place along this trajectory with the 



effective potential given by Vcs{4>) 
given by 

_ 1 / m^Mpi^ 

^ " 2 Vw4+ 1^202 



+ m?(fP'/2). Then, the slow-roll parameters are 



0. 



(109) 
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Figure 10: Contours of /nl (left panel), tnl (center panel) and (^nl (right panel) in the 
s-a^: plane for the chaotic inflation model with n = 6. The shaded region is excluded by 
WMAP3 data. 



For the parameters above, without the curvaton contribution, the number of e-folding, 
the spectral index and the tensor-to-scalar ratio are given by iVjnf = 55.3, = 0.975 
and r = 0.118. The running of the spectral index is negligible in this case. Although 
these values are consistent with the WMAP3 data, to see how the curvaton affects these 
predictions, contours of Ug and r are shown in Fig. [121 Furthermore, contours of non- 
linearity parameters are also shown in Fig. [131 In this case, large non-Gaussianity can also 
be generated. 

In the above choice of the parameters, we have red-tilted power spectrum in the absence 
of the curvaton. We can choose another set of parameters such that m <^ gv. In this case, 
we have t] ^ e and the spectrum becomes blue, i.e., Ug ^ 1 + Gt], in the absence of the 
curvaton. If the curvaton contribution is small enough to satisfy eM^^Q"^ <C 1, then adding 
the curvaton yields no effects on the parameters as in the case of the new inflation. On 
the other hand, if the curvaton contribution is large enough to satisfy eM^^Q"^ ^ 1, then 
the spectral index is given by - 1 ^ ~2e + 7]/{eM^^Ql). Because e is very small (e <^ 77) 
and f] is positive in this model, the spectrum is almost scale invariant or blue-tilted, and r 
is negligibly small. Hence, for m gv in which case the spectrum is blue with negligibly 
small scalar-to-tensor ratio, the model with such parameters is hardly allowed by WMAP3 
(see Fig. Hj) even if the curvaton contributes to the primordial fluctuations. 



4 Summary and Discussions 

In this paper, we studied primordial fluctuations in models where both the inflaton and the 
curvaton can contribute to the cosmic density fluctuations. By using the 6N formalism, 
we provided a systematic formulation to discuss the primordial curvature fluctuations and 
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Figure 11: Contours of in the s-a* plane for the new inflation model with n = 3 (left) 
and n = 4 (right). We take v = 10~^Mpi. In the absence of the curvaton contribution, 
Us = 0.920 and 0.944 for the case with n = 3,iVinf = 50.3 and n = 4, iVinf = 53.3 
respectively. 



their non-Gaussianity in this kind of scenario. We gave the expressions for the spectral 
index Ug, its running rirun, the tensor-to-scalar ratio r and the non-linearity parameters 
/nl, ttssl and (^nl for the mixed model. In general, the values of these parameters can be 
obtained by numerically calculating some set of equations, however, in some cases, they are 
given analytically. We presented such analytic formulae in some limiting cases. We found 
that very large non-Gaussianity (/nl ^ 1) can be generated even if the contribution 
of the curvaton fluctuations to the total ones is comparable to or smaller than that of 
the infiaton fluctuations. If the curvaton fluctuations are minor components, then tnl 
is enhanced by the factor of Ciif/Ccur ^ compared to /^l, where Cinf and Ccur are 
the primordial fluctuations from the infiaton and the curvaton respectively. In such a 
case, large non-Gaussian feature appears through the trispectrum rather than through 
the bispectrum. We also derived the consistency relation between the bispectrum and the 
trispectrum which holds in the limit of the large non-Gaussianity. If the measurement of 
the trispectrum becomes available in the future, this relation will provide a useful tool to 
distinguish the mixed model of the infiaton and the curvaton from other scenarios which 
also generate large non-Gaussianity. 

After we worked out the formulation for the calculation of inflationary parameters 
including the non-linearity parameters, we investigated the issue for several concrete mod- 
els of inflation: chaotic inflation, new inflation and hybrid inflation. By calculating the 
spectral index and the tensor-to-scalar ratio, we can compare those with cosmological ob- 
servations. We showed that some inflation models already excluded by the current data 
such as chaotic inflation with higher polynomials n = 4 and 6 can be liberated depending 
on the parameters which characterize the curvaton. By presenting contours of Ug and r in 
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Figure 12: Contours of (left panel) and r (right panel) in the s-a* plane for the hybrid 
inflation model with a quadratic potential. In the absence of the curvaton contribution, 
the spectral index and tensor-to-scalar ratio become as Ug = 0.975 and r = 0.118 for 
A^inf = 55.3. 



the s-cr* plane, we have explicitly shown in what cases such models can be relaxed to be 
allowed by the data. We have also presented contours of the non-linearity parameters in 
the s-(j* plane. In particular, it should be noticed that, the non-Gaussianity can be large 
due to the contribution from the curvaton without conflicting with the data. 

On the other hand, some inflation models such as the new inflation where the slow- 
roll parameter e is very small are scarcely affected by the curvaton. The effects of the 
curvaton on Ug and r appear as in the combination of eQ'^, thus the size of e is crucial 
when one considers how the curvaton affects the inflationary predictions. Furthermore, 
non-Gaussianity cannot be large in this model even if we introduce the curvaton because 
the e parameter is very small in this model. As shown in Eqs. P8l) ^ fl50|) . the curvaton 
contribution to the non-linearity parameters is always associated with e. Thus inflation 
models with very small e cannot generate large non-Gaussianity even with the curvaton. 

Although we have concentrated on adiabatic fluctuations thus far, isocurvature fluc- 
tuations can be produced in principle, depending on the thermal history of the universe. 
It can happen when dark matter or baryon was generated (and its number is conserved) 
before the decay of the curvaton or from the decay products of the curvaton [59l [601 EI] • 
In the case where dark matter and baryon were produced after the decay of the curvaton, 
no isocurvature perturbation is generated even in our mixed model. In this paper, we 
have considered this kind of case to avoid complexities coming from the contamination of 
isocurvature modes when we interpret observational signatures. However, here we com- 
ment on the cases where some isocurvature fluctuations are generated. First of all, since 
pure cold dark matter (CDM)/bayron isocurvature modes (regardless of its correlation 
with adiabatic mode) are severely constrained by observations, the pure curvaton scenario 
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Figure 13: Contours of /nl (left panel), tnl (center panel) and (^nl (right panel) in the 
s-a^: plane for the hybrid inflation model with a quadratic potential. 

with isocurvature fluctuations such as the above mentioned cases are disfavored from cur- 
rent observations. However, it should be noted here that the constraint can be relaxed 
when adiabatic fluctuations from the infiaton contribute to the total density fluctuations, 
which is the case for the mixed scenario. This is simply because the contribution from 
adiabatic mode decreases the fraction of the isocurvature contribution to the total fluc- 
tuations. Thus the mixed scenario may have interesting implications on scenarios with 
isocurvature fluctuations. For example, even in the case where dark matter was generated 
before the decay of the curvaton, the mixed scenario is still allowed if adiabatic fluctuations 
from the infiaton is bigger than that from the curvaton to some extent. Similarly, although 
a scenario where dark matter was generated from the decay products of the curvaton is 
severely constrained by the argument of isocurvature fluctuations, as in the case discussed 
above, the constraint is significantly relaxed in the mixed scenario if adiabatic fluctuation 
from the infiaton is much bigger than that from the curvaton. It should be noticed that 
large non-Gaussianity would not be produced when the contribution from the infiaton is 
large compared to that from the curvaton. On the other hand, when the contribution from 
the curvaton with isocurvature mode is large, such a scenario is severely constrained by 
the consideration of isocurvature mode, which corresponds to the case where the fraction 
of energy density of the curvaton to the total one becomes comparable to that of the 
infiaton. Remind that large non-Gaussianity can be generated when the fraction of energy 
density of the curvaton to the total one is very small. Thus non-Gaussianity would not be 
so large in this case too. 

Since a simple single- field inflation model can only generate very small non-Gaussianity, 
if the evidence of non-Gaussianity, which was recently reported in Ref . [3] , is established 
in the future, this mechanism of generating large non-Gaussianity will be very interesting. 
Although we concentrate on the mixed model of the infiaton and the curvaton in this paper, 
other mixed model may also be worth investigating. In particular, large non-Gaussianity 
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can also be generated in the modulated reheating [621 [631 EH [27] and preheating scenarios 
[651 [66l [67] so that the implications of such a mixed model on the inflationary parameters 
and non-Gaussianity are of great interest and will be presented elsewhere [68] . 
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Appendix 



A Value of a at the onset of the curvaton oscillations 

In the main text, we defined g{a^:) with the energy density of the curvaton at the onset of 
the oscillations as 

p.(iV^) = ^mlgia^y, (110) 

with g{a^) = aa^. We assume that the universe is radiation-dominated when the curvaton 
begins to oscillate, thus the equation of motion for a can be written as 

a + + mla = 0, (111) 

where we used H = 1/ (2t) during radiation-dominated epoch. The solution for this 
equation with the initial value (j(tm) = (y* is given by 

= '''|^^ff"* ^i/4M), (112) 

where J is the Bessel function. Using the asymptotic form of J1/4, the energy density 
for mt 3> 1 becomes 

Thus a can be given by 

a = ^r(5/4) ~ 1.45. (114) 

'71 



B Analytic expression of F(p) for ^ 1 

In this Appendix, we give a detailed description how we obtain an analytic expression of 
F{p) for the case with p <^ 1. We also assume that the value of s = Tcr/rricr is also very 
small. F{p) is given as 

F{P)= rdNe^-^P^. (115) 
Jo H{N) Pro 

To evaluate F{p), we have to follow the evolution of energy densities of the curvaton and 
radiation after the curvaton begins to oscillate for the case we are going to consider here. 
Thus we take the initial time = when H is equal to m„ and the curvaton field can be 



33 



considered to be dust. In this case, the background equations we need to follow are: 

^Pr-+V = ^Pa, (116) 

J^P. + 3P. = -^P., (117) 

^' = ^(Pr+P.), (118) 

where we use the number of e-folding as a time variable. We solve this set of equations 
with the initial conditions: 

Pr.(0)=p.o, Pa(0) = p.o = f^. H{0) = Ho = m^. (119) 

We assume that curvaton energy density is negligibly small throughout the evolution. 
Hence we treat Po-o as the expansion parameter of order p, where p is infinitesimal. Then 
we can expand p^, Pa,H with respect to p as 

Pr = + P« +Pf^ + ---. (120) 

P. = pi°^+P«+Pi^^ + ---, (121) 
H = + //(I) + + . . . ^ (122) 

where p^^\ p^a\ iJ^"^ = 0{p'^). We also expand F{p) with respect to p as 

F = + + + . . . . (123) 

Obviously, F^^^ vanishes. We are going to obtain the expression for F(p) by solving the 
above set of equations order by order in the following. 

• 0th order in p 

Since the energy density of the curvaton field is considered to be the order of 0{p), at the 
0-th order, it vanishes as plj^ = 0. Other quantities can be easily obtained. 

P^°^ = P.oe-^^, = Hoe-''^. (124) 

• 1st order in p 

From Eq. (11171) . the equation for p^a^ can be written as 

^pS^^+3pS^^ = -.e^^pS^). (125) 
By solving this equation, the energy density of the curvaton in the 1st order is given by 

(126) 



Pa ^ = PaO exp 



2 
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At the 1st order in p, F{p) can be given by 



dNe 



AN 



r 

J- cr l->cr 



(127) 



By substituting the expressions for p"a^ and if '•"^ into the above equation, we obtain 



oo p 



dNe " ' — ' „ exp 



-3N 



r 



(e^^-l) 



2Ha 



(128) 



The integral which appears in this equation can be done analytically as, 



dA^exp 



3N 



2m, 



(e- - 1) 



^(v^+v/^e^/2(l-Erf(v^))). (129) 



where Erf (x) is the error function and defined by 

Erf(a;) = ^ 



dte 



(130) 



By expanding the error function with s and taking the leading order, we can obtain F at 
the first order in p as 



(131) 



To evaluate the function F{p) at the second order, we need the expressions for pi^^ and 
H^^l For the energy density of radiation for the 1st order, from Eq. flll6p . we have the 
equation for pi^^ as 



dN^' +4Pr Pa ■ 

Using Eq. (11261) . the solution of this equation can be formally written as 



dN' exp(3iV'-^(e2^'-l) 



(132) 



(133) 



By plugging this expression into the Friedmann equation, the Hubble parameter can be 
given by, in the 1st order. 



1 



PcrO S. 

62 



//(o) 2 p,o 



s J dN' exp (SN' - '^e^^') + exp 



N - -e^^ 
2 



(134) 



• 2nd order in p 

Now we calculate the function F{p) for the second order in p. F^'^\p) can be written as 



dN e 



6N ro" 



-s dNe 



Pro Jo 
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Pro 



(135) 



To evaluate the RHS, we do not need pr^^ and iJ*^^-* because these do not appear in F{p) 
up to second order in p, but we need to have the expression for p^^\ From Eq. (ITTTl) . we 
have an equation for plj^ as 



_^ (2) , o (2) 



tt{1) 

Pa ^ H{0)P^ ■ 



(136) 



(2) 

Using the first order solutions obtained above, pa can be written as 



P. 



(2) _ ^ PcrO 



dN' e 



t: — p^oe exp -3iV - -e / 

/ PrO V Z / Jq 

/ dN" exp (3iV" - |e2^") + exp [n' - ^-e 

J 



(137) 



In principle, F^^^ can be obtained by plugging Eqs. (11331) . (I134p and (I137P into Eq. (11351) . 
the expression becomes lengthy and complicated, we consider each term in the RHS of 
Eq. (11351) in each. 

We denote the first term in Eq. (I135p as T\. Using the formulae for p^}^ and i/*^^-', T\ 
can be written as 



: / dN" exp (3iV" - ^e^^ 
JO V 2 



dW e 



I „2Ar' 



exp [n' - |e2^' 



(138) 



Although we need to perform a triple integral, we can make use the following formulaic 



exp (3iV' - |e-) . 0) - i^exp (iV - £e- 



4 



Erf 



(140) 



where we have neglected higher order terms in s because we are interested in the case 
s <C 1. Using this formula, the first term in Eq. (I138p can be written as Pi/s^ (omitting 
the factor appearing outside the integral and taking only leading order term in s) where 
Pi is given by 



8 ^ rfy y^e-y' dx (^-x'e""' + ^xErf (x)^ . (141) 



This equation can be verified by changing the variable as 



(139) 
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After some calculations, we find (3i = 1. The second term in Eq. fll38p can be calculated 
as 7r/(4s^). Hence Ti is given by 

Next we evaluate the second term in Eq. (11351) which is denoted as T2. By using the 
expressions for pi-^-* and H^^\ T2 can be written as 



s J dN' exp (3iV' - |e2^') + exp ( 



(143) 



Similar calculations to the case with Ti yields 



1 . s / \ 2 



where we have again dropped higher order terms in s. Thus we obtain the expression for 

as 

=T,-T2 = (145) 
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